RHESSI and SOHO/CDS Observations of Explosive 
Chromospheric Evaporation 
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ABSTRACT 

Simultaneous observations of explosive chromospheric evaporation are presented using data 
from the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) and the Coronal 
Diagnostic Spectrometer (CDS) onboard SOHO. For the first time, co-spatial imaging and spec- 
troscopy have been used to observe explosive evaporation within a hard X-ray emitting region. 
RHESSIX-r&y images and spectra were used to determine the flux of non-thermal electrons accel- 
erated during the impulsive phase of an M2.2 flare. Assuming a thick-target model, the injected 
electron spectrum was found to have a spectral index of ~7.3, a low energy cut-off of ~20 keV, and 
a resulting flux of >4xl0 10 ergs cm~ 2 s _1 . The dynamic response of the atmosphere was deter- 
mined using CDS spectra, finding a mean upflow velocity of 230±38 km s" 1 in Fe XIX (592.23 A), 
and associated downflows of 36±16 km s _1 and 43±22 km s _1 at chromospheric and transition 
region temperatures, respectively, relative to an averaged quiet-Sun spectra. The errors repre- 
sent a lcr dispersion. The properties of the accelerated electron spectrum and the corresponding 
evaporative velocities were found to be consistent with the predictions of theory. 

Subject headings: Sun: atmospheric motions - Sun: flares - Sun: UV radiation - Sun: X-rays, 7 rays 



1. INTRODUCTION 

Current solar flare models (Antiochos & Stur- 
rock 1978; Fisher, Canfield, & McClymont 1984, 
1985a,b,c; Mariska, Emslie, & Li 1989) predict 
two types of chromospheric evaporation processes. 
"Gentle" evaporation occurs when the chromo- 
sphere is heated either directly by non-thermal 
electrons, or indirectly by thermal conduction. 
The chromospheric plasma subsequently loses en- 
ergy via a combination of radiation and low- 
velocity hydrodynamic expansion. "Explosive" 
evaporation takes place when the chromosphere is 
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unable to radiate energy at a sufficent rate and 
consequently expands at high velocities into the 
overlying flare loops. The overpressure of evapo- 
rated material also drives low-velocity downward 
motions into the underlying chromosphere, in a 
process known as chromospheric condensation. 

From a theoretical perspective, Fisher, Can- 
field, & McClymont (1985a) investigated the re- 
lationship between the flux of non-thermal elec- 
trons (F) and the velocity response of the atmo- 
sphere for the two classes of evaporation. For gen- 
tle evaporation, non-thermal electron fluxes of < 
10 10 ergs cm~ 2 s _1 were found to produce upflow 
velocities of tens of kilometres per second. In con- 
trast, explosive evaporation was found to be as- 
sociated with higher non-thermal electron fluxes 
(F > 3 x 10 10 ergs cm~ 2 s _1 ) which drive both 
upflows of hot material at velocities of several 
hundred kilometres per second and downflows of 
cooler material at tens of kilometres per second. 

Observationally, previous studies have identi- 
fied blue-shifted Soft X-Ray and EUV lines indica- 
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Fig. 1 — CDS images in the He I, O V, and Fe XIX 
emission lines observed during the implusive phase 
of the flare, with the corresponding EIT 195 A im- 
age. RHESSI 12-25 keV (dashed) and 25-60 keV 
(solid) contours are overlayed, each drawn at 10% 
of the peak intensity. 



tive of chromospheric evaporation. Using the Bent 
Crystal Spectrometer onboard the Solar Maxi- 
mum Mission, Antonucci & Dennis (1983) and 
Zarro & Lemen (1988) reported upflow velocities 
of 400 km s _1 and 350 km s -1 , respectively, in 
Ca XIX lines (3.1-3.2 A). More recently, Cza- 
ykowska et al. (1999), Teriaca et al. (2003), and 
Del Zanna et al. (2005) observed velocities of 140- 
200 km s" 1 in Fe XIX (592.23 A), using the Coro- 
nal Diagnostic Spectrometer (CDS; Harrison et al. 
1995) onboard the Solar and Heliospheric Obser- 
vatory (SOHO). Simultaneous upflows and down- 
flows during a Hard X-Ray (HXR) burst indica- 
tive of explosive evaporation have been observed 
using CDS and Yohkoh/Haxd X-Ray Telescope by 
Brosius & Phillips (2004). While these studies 
provided a measurement of the dynamic response 
of the flaring chromosphere, they were unable to 
provide a measurement of the flux of electrons re- 
sponsible for driving such motions, nor the spatial 
relationship between the two. 

In this Letter, simultaneous Reuven Ramaty 
High Energy Solar Spectroscopic Imager (RHESSI; 



Lin et al. 2002) and CDS observations are com- 
bined for the first time to investigate the relation- 
ship between the non-thermal electron flux and the 
response of the solar atmosphere. In Section 2 the 
analysis techniques employed are described, while 
the results are presented in Section 3. Our Con- 
clusions are then given in Section 4. 

2. OBSERVATIONS AND DATA ANAL- 
YSIS 

This study focuses on a GOES M2.2 flare, which 
began at 12:44 UT on 2003 June 10. The event was 
selected from a sample of approximately 50 flares 
jointly observed by RHESSI and CDS. The limited 
field of view, cadence, and operating schedule of 
CDS, coupled with RHESSI nighttime and South 
Atlantic Anomaly passes, make simultaneous ob- 
servations by the two instruments quite rare. 

2.1. The Coronal Diagnostic Spectrome- 
ter (CDS) 

The CDS observations reported here were ob- 
tained with the FLARE_AR observing sequence. 
FLARE_AR contains five <4 A wide spectral win- 
dows centered on He I (584.33 A; log T = 4.5), O V 
(629.73 A; log T = 5.4), Mg X (624.94 A; log T = 
6.1), Fe XVI (360.76 A; log T = 6.4), and Fe XIX 
(592.23 A; log T = 6.9). Each raster consists of 
45 slit positions, each ~15 seconds long, result- 
ing in an effective cadence of —11 minutes. The 
slit itself is 4"xl80" yielding a -180" x 180" field 
of view. A zoomed-in region of the He I, O V, 
and Fe XIX rasters from the impulsive phase of 
the flare is given in Figure 1. Also shown is the 
Extreme ultraviolet Imaging Telescope (EIT; De- 
laboudiniere et al. 1995) 195 A passband image ob- 
tained at 12:48 UT. A series of subsequent EIT im- 
ages makes it clear that the He I and O V bright - 
enings come from a flare ribbon rather than two 
distinct footpoints as Figure 1 may suggest. 

The spectrum from each CDS pixel was fit- 
ted with a broadened Gaussian profile (Thomp- 
son 1999), for each of the five spectral windows. 
Velocity maps were created by measuring Doppler 
shifts relative to quiet-Sun spectra, which were as- 
sumed to be emitted by stationary plasma. Pre- 
liminary fits to the Fe XIX line during the im- 
pulsive phase of the flare revealed an asymmetric 
broadening beyond the instrumental resolution of 
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Fig. 2. — Sample of spectra from the spectral window centered on Fe XIX (592.23 A). The vertical short- 
dashed line shows the rest wavelength of the Fe XIX line at 592.32 A while the long-dashed line shows the 
rest wavelength of the Fe XII line at 592.62 A. Panel 'a' was obtained from a quite Sun region, panel 'b' 
shows a stationary Fe XIX line from a post-implusive phase flare kernel, while panel 'c' shows an Fe XIX line 
from the flare ribbon during the implusivc phase. The dotted line indicates the stationary component while 
the triple-dot-dashed line indicates the blue-shifted component. 



CDS. The strongest blue asymmetries were found 
within the flare ribbon during the impulsive phase. 
Outside this area, and after the impulsive phase, 
the Fe XIX line was observed to have a width 
comparible to the instrumental width. Figure 2 
shows a sample of spectra taken from the spec- 
tral window centered on the Fe XIX (592.23 A) 
emission line. Panel 'a' shows a spectrum from 
a quiet Sun area in which no Fe XIX emission 
was visible. Instead, a weak emission line was 
observed at 592.6 A which Del Zanna & Mason 
(2005) have identified as Fe XII. Panel 'b' shows 
a stationary Fe XIX emission line extracted from 
a bright region, but after the impulsive phase at 
^12:50 UT when no significant flows are expected. 
An emission line with a strong blue asymmetry is 
shown in panel 'c'. This was extracted from the 
flare ribbon during the impulsive phase. The best 
fit to this line was consistent with stationary and 
blue-shifted components, both with widths com- 
parible to the instrumental resolution. As Fe XIX 
is not observed in quiet-Sun spectra, and follow- 
ing from Teriaca et al. (2003), the Doppler veloc- 
ity was measured as the shift between these two 
components. A hcliographic correction was also 
applied, due to the longitude of the observations 
and assuming purely radial flows, 



2.2. The Reuven Ramaty High Energy 
Spectroscopic Imager (RHESSI) 

RHESSI is an imaging spectrometer capable of 
observing X- and 7-ray emission over a wide range 
of energies (~3 keV-17 MeV). During this event 
the thin attenuators on RHESSI were in place 
thus limiting the energy range to >6 kcV. Flare 
emission was not observed above ~60 keV. The 
flare lightcurves are shown in the top panel of Fig- 
ure 3. Both the RHESSI images and spectra were 
obtained over a 64 second period from 12:47:34- 
12:48:38 UT to coincide with the timerange over 
which CDS observed blue asymmetries in the 
Fe XIX line. This time interval lies within the 
impulsive 25-60 keV HXR burst and is indicated 
by two vertical dotted lines in the top panel of 
Figure 3. RHESSI images in two energy bands 
(12-25 and 25-60 keV) were reconstructed using 
the Pixon algorithm (Hurford et al. 2002). Con- 
tours at 10% of the peak intensity in each band 
are overlayed on each EUV image in Figure 1. 

The RHESSI spectrum was fitted assuming an 
isothermal distribution at low energies, and thick- 
target emission at higher energies (bottom panel 
of Figure 3). A thick-target model was chosen over 
a thin-target model as it is believed that the den- 
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Fig. 3. — Top panel: RHESSI lightcurvcs from 
the 6-12, 12-25, and 25-60 keV bands. The dot- 
ted vertical lines indicate the time interval over 
which images and spectra were obtained to corre- 
spond to time when significant upflows were ob- 
served using CDS. Bottom panel: Portion of the 
RHESSI spectrum integrated over the timcrange 
given above. The energy range 6-60 keV (ver- 
tical dot-dash lines) was fitted with an isother- 
mal component (dotted curve) and a thick-target 
bremsstrahlung component (dashed curve). 



sity of the flare loop is insufficent to thcrmalisc 
the electrons as they propagate to the chromo- 
sphere. The thick-target model is used in the vast 
majority of cases (e.g. Holman 2003, Veronig et 
al. 2005). Furthermore, in this flare the HXR 
source is clearly aligned with the He I ribbon 
as seen by CDS, which implies that the acceler- 
ated electrons are losing their energy in the dense 
chromosphere rather than in the coronal loops. 
The total power of non-thermal electrons above 
the low energy cut-off (e c ) was calculated from 
P(e > e c ) = f e (e)de ergs s _1 , where f e (e) ~ 



e~ s electrons keV" 1 s _1 is the thick-target elec- 
tron injection spectrum and S is the associated 
spectral index (Brown 1971). Because of the steep- 
ness of the RHESSI spectrum at high energies, 
the non-thermal flux is quite sensitive to the value 
of the low energy cut-off. In order to put a con- 
straint on this value, the temperature of the ther- 
mal component was obtained by another indepen- 
dent method, i.e. the equivalent width of the Fe 
line complex at 6.7 keV (Phillips 2004). The value 
of the equivalent width of this line, which is quite 
sensitive to the temperature, was used to estimate 
the temperature of the thermal component. Hav- 
ing fixed this value, the entire RHESSI spectrum 
was fitted using a least-squares fit. 

3. RESULTS 

The thick-target model fitted to the RHESSI 
spectrum in Figure 3 was consistent with an elec- 
tron distribution having e c ~20 keV and 6 <~7.3. 
The break energy of 20 keV is consistent with 
earlier works (e.g. Holman 2003, Sui, Holman, 
& Dennis 2005). The total power in non-thermal 
electrons was therefore lxlO 29 ergs s _1 . Explor- 
ing the possible range of values for the break en- 
ergy for this flare would yield an electron power 
value of 4xl0 29 ergs s _1 for e c = 15.0 keV 
while e c = 25.0 keV would give a power value of 
6xl0 28 ergs s _1 . However, either of these break 
energies would give a worse \ 2 value than ob- 
tained from the original fit. By comparison, the 
total thermal power for the same time interval was 
found to be 1.2xl0 28 ergs s _1 . 

Using the reconstructed 25-60 keV image, an 
upper-limit to source size was calculated to be 
2.3 xlO 18 cm 2 . This was found by summing over 
all pixels with counts greater than 10% of the peak 
value. This threshold was chosen to eliminate 
sources outside of the main HXR-emitting region, 
which were assumed to be unreal; the source area 
was not found to be highly sensitive to this value. 
For example, a threshold of 5% yielded an area of 
3.2 xlO 18 cm 2 and 15% yielded 2xl0 18 cm 2 . This 
area was also confirmed using the Fourier mod- 
ulation profiles from each of RHESSIs nine de- 
tectors, which are sensitive to spatial scales from 
2.2"to 183". Assuming a filling factor of unity, the 
resulting non-thermal electron flux was calculated 
to be >4xl0 10 ergs cm~ 2 s _1 . 
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Fig. 4. — Velocity maps in He I and Fe XIX. 
Downflows are indicated by red pixels while up- 
flows are indicated by blue pixels. The vertical 
dashed lines correspond to the times indicated by 
the vertical dashed lines in Figure 3 and the ar- 
row denotes the direction in which the CDS slit 
moves. Black regions in the Fe XIX map represent 
pixels where no significant Fe XIX emission was 
observed. RHESSI 25-60 keV contours at 10% 
and 40% of the peak intensity are overlayed. 



Figure 4 shows velocity maps in the He I and 
Fe XIX lines. The He I map shows consistent 
downflows of 20-50 km s" 1 until the slit leaves 
the flaring region at —12:50 UT. A velocity map 
in O V showed a similar trend. However, the 
Fe XIX map shows strong upflow velocities of 190- 
280 km s" 1 during the HXR peak, indicated by 
'Ti' and 'T2' on Figures 3 and 4. No significant 
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Fig. 5. — Plasma velocity as a function of temper- 
ature for the five lines observed using CDS. Pos- 
itive velocities indicate downflows, while negative 
values indicate upflows. 



upflows were evident once the HXRs begin to di- 
minish from time T2 onwards. 

By identifying each pixel in the Fe XIX map 
that required a two-component fit to the line pro- 
file, between times Ti and T2, the velocity was 
measured for the corresponding pixel in each of 
the five CDS rasters. Figure 5 shows the mean ve- 
locity as a function of temperature for each line us- 
ing the methods described in Section 2. The error 
bars represent a la dispersion. At chromosphcric 
and transition region temperatures, plasma ve- 
locities show red-shifts of 36±16 km s" 1 and 
43±22 km s , respectively, while the blue-shift 
observed in the 8 MK Fe XIX line corresponds 
to a velocity of 230±38 km s _1 . No significant 
flows were observed in the Mg X and Fe XVI lines. 
The combination of high- velocity upflows and low- 
velocity downflows, together with a non-thermal 
electron flux of >4xl0 10 ergs cm" 2 s" 1 provides 
clear evidence for explosive chromosphcric evapo- 
ration. 

4. DISCUSSION AND CONCLUSIONS 

For the first time, co-spatial and co-temporal 
HXR and EUV observations of chromospheric 
evaporation are presented using RHESSI and 
SOHO/ CDS. High upflows velocities (-230 km s" 1 ) 
were clearly observed in high-temperature Fe XIX 
emission during the impulsive phase of an M2.2 
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flare, while much lower downflow velocities (^40 km s 
were observed in the cooler He I and O V 
lines. The value of the non-thermal electron flux 
(>4xl0 10 ergs cm~ 2 s _1 ) and the resulting veloc- 
ity response are indicative of an explosive evapo- 
ration process occuring during this flare, as laid 
out in Fisher, Canfield, & McClymont (1985a) 
and Mariska, Emslie, & Li (1989). 

The combination of HXR and EUV observa- 
tions presented in this Letter have enabled us 
to obtain a greater understanding of the char- 
acteristics of chromospheric evaporation, a fun- 
damental process in solar flares. We have pre- 
sented the first detection of explosive mass mo- 
tions within HXR footpoints, and determined the 
flux of non-thermal electrons responsible for driv- 
ing such flows. 
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